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43 The MiniBooNE and SciBooNE collaborations report the results of a joint search for short base- 

44 line disappearance of at Fermilab's Booster Neutrino Beamline. The MiniBooNE Cherenkov 

45 detector and the SciBooNE tracking detector observe antineutrinos from the same beam, therefore 

46 the combined analysis of their datasets serves to partially constrain some of the flux and cross sec- 

47 tion uncertainties. Uncertainties in the background were constrained by neutrino flux and cross 

48 section measurements performed in both detectors. A likelihood ratio method was used to set a 

49 90% confidence level upper limit on P M disappearance that dramatically improves upon prior limits 
so in the Am 2 =0. 1-100 eV 2 region. 
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I. INTRODUCTION 33 

34 

Recently there has been increasing evidence for neu- 35 
trino oscillations in the Am 2 sa 1 eV 2 region. Such os- 36 
dilations could imply the existence of new sterile neu- 37 
trino species that do not participate in standard weak 38 
interactions but mix with the known neutrino flavors 39 
through additional mass eigenstates. The LSND [l[ ex- 40 
periment observed an excess of j> e -like events in a i/„ 41 
beam. MiniBooNE has observed an excess of z/ R -like 42 



and zA,-like events, 



in a Vfj, beam and 



beam, respec- 



11 tively. Additional evidence for short-baseline anomalies 44 

12 with L/E w 1, where L is the neutrino path length in km 45 

13 and E the neutrino energy in GeV, includes the deficit of 46 

14 events observed in reactor antineutrino Q and radioac- 47 

15 tivc source neutrino @ experiments. If these anomalies 48 

16 are due to neutrino oscillations in the Am 2 ~ 1 eV 2 49 

17 range, then sterile neutrinos may provide the additional 50 
is mass eigenstate(s). Observation of (v^) disappearance 51 

19 in conjunction with v e (D e ) appearance in this Am 2 range 52 

20 would be a smoking-gun for the presence of these sterile 53 

21 neutrinos. Alternatively, constraining disappear- 54 

22 ance can, along with global v e (v e ) disappearance data, 55 

23 constrain the oscillation interpretation of the v e (v e ) ap- 56 

24 pearance signals in LSND and MiniBooNE Q. 57 

25 Searches for and disappearance in MiniBooNE 58 

26 were performed in 2009 ||. No evidence for disappear- 

27 ance was found. The search for disappearance was 

28 recently repeated in MiniBooNE with the inclusion of 59 

29 data from the SciBooNE detector in a joint analysis @. 60 

30 Once again, the results were consistent with no dis- 

31 appearance. The analysis presented here is an improved ei 

32 search for disappearance using data from MiniBooNE 6 2 
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and SciBooNE taken while the Booster Neutrino Beam- 
line (BNB) operated in antineutrino mode. 

The Monte Carlo (MC) predictions for both Mini- 
BooNE and SciBooNE were updated to account for re- 
cent neutrino flux and cross section measurements made 
with both experiments. The data from both detectors 
were then simultaneously fit to a simple two-antineutrino 
oscillation model. Improved constraints on MC predic- 
tions, the inclusion of SciBooNE data, and a MiniBooNE 
antineutrino data set nearly three times larger than what 
was available for the original disappearance analysis, 
have allowed a 90% confidence level upper limit to be set 
that dramatically improves upon prior limits in the Am 2 
= 0.1-100 eV 2 region, pushing down into the region of 
parameter space of interest to sterile neutrino models. 

This paper is organized as follows. Section [TT1 describes 
the BNB and the MiniBooNE and SciBooNE detectors. 
Then, the simulation of neutrino interactions with nu- 
clei and subsequent detector responses arc described in 
Section Mil The event selection and reconstruction for 
both detectors are described in Section HVl The param- 
eters for the MC tuning and its systematic uncertainties 
are given in Section [Vj Section IVTl describes the analysis 
methodology. The results of the analysis are presented 
in Section IVII1 and the final conclusions are given in Sec- 
tion [vm] 



II. BEAMLINE AND EXPERIMENTAL 
APPARATUS 

MiniBooNE and SciBooNE both use the BNB at Fer- 
milab in Batavia, Illinois. The 8 GeV kinetic energy pro- 
tons from the booster accelerator strike a 1.7 interaction 
length beryllium target, which is located inside a focus- 
ing horn. The horn is pulsed in time with the beam to 
produce a toroidal magnetic field that, depending on the 
polarity setting, will either focus n~/K~ and defocus 
ir Jr /K + or vice- versa. These mesons then pass through 
a 60 cm long collimator and decay in flight along a 50 
m long tunnel. A schematic view of the BNB from the 
beryllium target to both detectors is shown in Figure [T] 

The resulting neutrino beam will have an enhanced 
flux of either muon neutrinos (neutrino mode) or muon 
antineutrinos (antineutrino mode). In antineutrino mode 
beam running, the flux of antineutrinos in the beam will 
be referred to as the right-sign (RS) flux and the flux of 
neutrinos in the beam will be referred to as the wrong- 
sign (WS) flux. These two designations are used because 
antineutrinos are the signal in this analysis and neutrinos 
are an intrinsic background. Figure [2] shows the neutrino 
and antineutrino flux prediction in antineutrino mode 
at both the MiniBooNE and SciBooNE detectors. De- 
tails on the beamline and flux predictions are given in 
Ref. [HI. 

The MiniBooNE detector (TTI is located 541 m down- 
stream of the antineutrino production target and consists 
of a spherical 12.2 m diameter tank containing 800 tons of 
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(a) MiniBooNE antincutrino mode flux 
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(b) SciBooNE antincutrino mode flux 



FIG. 2. The neutrino and antineutrino flux prediction as a function of true neutrino(antineutrino) energy, in antineutrino mode 
at the MiniBooNE and SciBooNE detectors. The flux is represented by the solid line, the flux is represented by the 
dashed line, the v e flux is represented by the dot-dashed line, and the v e flux is represented by the dotted line. 



1 mineral oil (CH2), beneath at least 3 m of earth overbur- 

2 den. The fiducial volume is a sphere 10 m in diameter, 

3 with a fiducial mass of 450 tons. The detector is instru- 

4 mented with 1280 8 inch photomultiplier tubes (PMTs) 

5 in the active region, and 240 8 inch PMTs in an outer, 
e veto region. Events are reconstructed based on timing 
7 and charge information mostly from Cherenkov radia- 
s tion. A schematic of the MiniBooNE detector is shown 

9 in Figure [3J 

10 The SciBooNE detector [l2| is located 100 m down- 

11 stream of the target. SciBooNE is a discrete tracking 

12 detector comprised of three sub-detectors (in order from 

13 upstream to downstream): a fully active and finely seg- 
« mented scintillator tracker (SciBar), an electromagnetic 

15 calorimeter (EC), and a muon range detector (MRD). 

16 The SciBar sub-detector [l3[ consists of 14336 extruded 

17 polystyrene (CsHg) strips arranged vertically and hori- 

18 zontally to construct a 3 x 3 x 1.7m 3 volume. Each scin- 

19 tillator strip is read out by a wavelength shifting fiber 

20 (WLS) attached to a 64-channel multi-anode PMT (MA- 




FIG. 3. Schematic view of the MiniBooNE detector. 
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PMT). The 15 ton SciBar sub-detector (10.6 ton fiducial 36 
volume) provides the primary interaction target. The 
EC sub-detector is a two plane (vertical and horizon- 
tal) "spaghetti" -type calorimeter; 64 modules made of 37 
1 mm scintillating fibers embedded in lead foil are bun- 3S 
died and read out at both ends by PMTs. The MRD 
sub-detector, designed to measure muon momentum, is 

. . 39 

made from 12 iron plates, each 5 cm thick, sandwiched 

. . ... 40 

between 13 alternating horizontal and vertical scintilla- 

. 41 

tor planes of thickness 6 mm that are read out via 362 

42 

individual 2 inch PMTs. A schematic of the SciBooNE 
detector is shown in Figure 0] 



SciBar 



MRD 




2m 



FIG. 4. Schematic view of the SciBooNE detector. 



III. MONTE CARLO SIMULATION 67 

68 

Simulation of the neutrino and antincutrino flux, neu- 69 
trino and antincutrino interactions in the detector, and 70 
detector response has been discussed in detail in previous 71 
publications for MiniBooNE [II [01 and SciBooNE [Hj]. 72 
Calculation of the neutrino and antineutrino flux at the 73 
detector is done with a GEANT4-based model [l?} that 74 
is constrained by external measurements [ltj EH and ac- 75 
counts for proton transport to the target, p-Be interac- 76 
tions in the target, meson production, focusing by the 77 
magnetic horn, meson propagation and decay, and neu- 78 
trino and antincutrino propagation to the detectors. 79 

Neutrino and antineutrino interactions in both detec- so 
tors arc simulated using the NUANCE [T|| event gener- ai 
ator. Bound nucleons are described by the Rclativistic 82 
Fermi Gas (RFG) model 0. The MiniBooNE detector a 3 
response is simulated using GEANT3 [2l| , which takes 84 
the final-state particles emerging from a nucleus and as 
propagates them through the detector. The GEANT3 86 
code was modified to include a custom model for light 87 
propagation in the detector [22| and to use GCALOR [23| as 
for pion absorption and charge exchange in the detector 89 
medium. SciBooNE uses GEANT4 [24| to simulate the 90 



interactions of hadronic particles with detector materials. 



IV. EVENT SELECTION AND 
RECONSTRUCTION 

MiniBooNE data from a total of 1.01 x 10 21 protons on 
target (POT) operation in antineutrino mode, from July 
2006 up through April 2012, are included in the analysis. 
Data from SciBooNE antincutrino mode operation from 
June 2007 through August 2008 are included, comprising 
a total of 1.53 x 10 20 POT for the SciBooNE contribution. 

MiniBooNE event selection and reconstruction is es- 
sentially identical to that used for a previous neutrino 
mode cross section measurement [14j. Events with 
only a single fj, + in the detector are selected. Event se- 
lection cuts are based on the beam timing, fiducial vol- 
ume, observation of two correlated events (the muon and 
its decay electron), and the likelihood of the fit to the 
muon hypothesis. These cuts are designed to reject in- 
coming particles (i.e. muons from cosmic rays or from 
neutrino and antincutrino interactions in the surround- 
ing material), ensure that the event is contained within 
the detector, and ensure correct event classification as 
well as accurate muon energy estimation. The capture 
of fx" resulting from initial CCQE interaction events 
is simulated in the MC and these specific events are not 
selected. In antineutrino mode, a sizeable fraction of the 
events (roughly 20%) are due to interactions. Mini- 
BooNE cannot distinguish between and P M events on 
an event- by- event basis, so /i~s from interactions are 
an irreducible background. 

For SciBooNE, the event selection and reconstruc- 
tion is nearly identical to the previous inclusive charged 
current measurement (l6| . Two-dimensional SciBar 
tracks are reconstructed using a cellular automaton al- 
gorithm [25| from SciBar hits. Three-dimensional SciBar 
tracks arc then reconstructed based on the timing and 
endpoint positions of the two-dimensional SciBar tracks. 
Two-dimensional tracks in the MRD are independently 
reconstructed using hits in the MRD that are clustered 
within a 50 ns timing window. Three-dimensional tracks 
in the MRD are reconstructed by matching the timing of 
the two-dimensional projections. If the downstream edge 
of a SciBar track lies in the last two layers of SciBar, a 
search for a matching track or hits in the MRD is per- 
formed. The upstream edge of the MRD track is required 
to be on cither one of the first two layers of the MRD, 
and to be within 30 cm of the projected entry point of the 
SciBar track into the MRD (a more detailed des crip tion 
of the track reconstruction can be found in Ref. jl2|). 

To select /i + events, the highest momentum track per 
event in the beam on-time window is required to have 
> 0.25 GeV/c to reduce the number of neutral current 
(NC) events. The energy loss of the track in SciBar must 
be consistent with a muon hypothesis, and must origi- 
nate within the 10.6 ton SciBar fiducial volume. These 
muon candidate tracks are further categorized as SciBar- 
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1 stopped or MRD-stoppcd. SciBar-stoppcd events have 55 

2 the downstream endpoint of the muon candidate track 56 

3 contained in the SciBar fiducial volume. MRD-stopped 57 

4 events have the muon candidate track being a SciBar 58 

5 track matched to MRD hits or to an MRD track with a 59 
e downstream endpoint that does not exit the back or sides eo 

7 of the MRD. Both SciBar-stopped and MRD-stopped 

8 events are used in the analysis. SciBooNE has no over- 

9 burden so cosmic backgrounds must be subtracted. For 

10 cosmic background estimation, the same muon selection 
u criteria are applied to a beam-off time window that is five 

12 times longer than the beam-on window. This event rate 

13 is scaled and subtracted from the beam-on data. 

14 The selected events include and interactions on 

15 carbon and hydrogen in the detectors. The reconstructed 

16 antincutrino energy is based on the assumption that the 

17 interaction is always a charged current quasielastic 
is (CCQE) interaction with a proton at rest in carbon: + 

19 p — > fi + n. Hence, it is a function of the measured energy 

20 and direction of the outgoing muon. The equation for 

64 

21 reconstructed energy is: 



2{M p -E b -E^ + P^cob6 ix ) ' 68 

(1)69 

22 where M n and M p are the mass of the neutron and pro- 70 

23 ton, Mfj_, E^, Pp, and 6^ are the mass, energy, momen- 71 

24 turn, and direction of the outgoing muon, and Eg is the 72 

25 binding energy (30 MeV for protons in carbon). Equa- 73 

26 tion [T] is applied to all events in data and MC, even 74 

27 though a sizeable fraction of the events are not CCQE 75 

28 (i.e. charged current single tt (CCItt) or multi-nucleon 76 

29 knockout events misidentified as CCQE events). Mini- 77 

30 B00NE has an estimated resolution for reconstructed en- 78 

31 ergy of 8.3% for CCQE events and 13.9% for all events. 79 

32 SciBooNE has a reconstructed energy resolution of 9.6% 

33 for CCQE events and 24.6% for all events. 

34 MiniBooNE and SciBooNE data and MC are put in so 

35 2Tbin histograms of E® E . The binning goes from 300 

36 MeV to 1.9 GeV, with individual bin widths as follows: ai 

37 bin 1, 100 MeV; bins 2-19, 66.7 MeV; bin 20, 100 MeV; a2 

38 bin 21, 200 MeV. The first and last two bins are wider 83 

39 to ensure adequate event statistics in data and MC. S4 

40 Figure [5] shows the predicted event distributions in as 

41 MiniBooNE for reconstructed antineutrino and neutrino m 

42 energy, for events on hydrogen and carbon nuclei. Fig- 8 7 

43 urc [5] shows the predictions for MiniBooNE's recon- 8 8 

44 structed antineutrino and neutrino energy distributions 89 

45 by interaction type: CCQE, CCItt, and all other inter- 90 

46 action types (CC multi-pion and NC). Table fl] shows the 91 

47 MC predictions for the selected MiniBooNE events by 92 

48 neutrino and interaction type. v e and v e contamination 93 

49 is negligible. 94 

50 The following plots show several properties of the se- 95 

51 lected SciBooNE events, as predicted by simulation. Fig- 96 

52 ure [7] shows the reconstructed antincutrino and neutrino 97 

53 energy distributions for events on hydrogen and carbon 98 

54 nuclei. Figure [5] shows the reconstructed antincutrino 99 



and neutrino energy distribution by interaction type: 
CCQE, CCItt, and other (CC multi-pion and NC). Ta- 
bleUshows the MC predictions for the selected SciBooNE 
events by neutrino and interaction type. The data set is 
estimated to contain an additional 811 events from cos- 
mic ray unions. 



TABLE I. MC predictions for the number of selected events 
by neutrino and interaction type in both MiniBooNE and 
SciBooNE. 





MiniBooNE 


SciBooNE 


interaction type 


v events 


v events 


v events 


v events 


CCQE 


37428 


9955 


4619 


1359 


CCItt 


8961 


2593 


1735 


1006 


CC mufti- tt or NC 


2364 


460 


959 


610 



The difference in shape between the SciBooNE RS 
(WS) and MiniBooNE RS (WS) energy distributions 
is mainly due to different event selection criteria be- 
tween MiniBooNE and SciBooNE. MiniBooNE selects 
for CCQE interaction events and SciBooNE selects for 
all CC interaction events so the SciBooNE sample has a 
larger percentage of non-CCQE interaction events. Since 
the antineutrino energy reconstruction is based on a 
CCQE interaction assumption, there are more SciBooNE 
events with a larger discrepancy between true antineu- 
trino energy and reconstructed antineutrino energy than 
in MiniBooNE, leading to shape differences. Differences 
in selection efficiency, antineutrino flux at the detector 
locations, and background rejection between MiniBooNE 
and SciBooNE also contribute to the shape differences. 

Figure |9] shows the distribution of the combined an- 
tincutrino and neutrino propagation distances, from pro- 
duction in the decay tunnel to interaction in SciBooNE 
or MiniBooNE. 



V. SYSTEMATIC UNCERTAINTIES 

Beam and cross section uncertainties are calculated 
for both MiniBooNE and SciBooNE using the multisim 
method [26j . In this procedure, groups of correlated sim- 
ulation parameters associated with beam production and 
cross section modeling uncertainties are sampled accord- 
ing to their covariance matrices. The parameters for each 
source of uncertainty (t^, K + , etc.) are sampled 1000 
times to obtain sufficient statistics. Each MC event in 
MiniBooNE and SciBooNE is reweighted based on these 
varied parameters forming 1000 new MC predictions of 
the E® E distribution in both detectors. Covariance ma- 
trices, in bins of E® E , are then computed for each source 
of uncertainty by comparing these 1000 new MC predic- 
tions to the default MC prediction. The procedure takes 
care of the correlation of beam production and cross sec- 
tion uncertainties between MiniBooNE and SciBooNE. 
Cross section and nuclear model uncertainties for v and 
v events are treated as uncorrelated due to the poor un- 
derstanding of differences between v and v interactions 
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(a) MiniBooNE RS Events 



(b) MiniBooNE WS Events 



FIG. 5. Reconstructed antineutrino and neutrino energy (E® E ) distributions for selected RS and WS MiniBooNE events on 
different target types (hydrogen or carbon) from MiniBooNE MC. Total events are represented by the solid line, events with 
interaction on carbon are represented by the dashed line, and events with interaction on hydrogen are represented by the 
dot-dashed line. 
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(a) MiniBooNE RS Events 



(b) MiniBooNE WS Events 



FIG. 6. Reconstructed antineutrino and neutrino energy (E® E ) distributions for selected RS and WS MiniBooNE events for 
different interaction types (CCQE, CCI7T, other) from MiniBooNE MC. Total events are represented by the solid line, CCQE 
interaction events are represented by the dashed line, CCI7T interaction events are represented by the dot-dashed line, and all 
other interaction (CC multi-7r or NC) events are represented by the short-dashed line. 



in nuclear modeling. Some detector specific uncertain- io 
ties are calculated using the unisim method [26| where u 
uncorrelated detector specific uncertainties are varied up 12 
or down by one standard deviation. 13 



A. Beam Uncertainties 

16 
17 

Uncertainties in the delivery of the primary proton is 
beam to the beryllium target, the primary beam optics, 19 
secondary hadron production in proton-beryllium inter- 20 
actions, hadronic interactions in the target and horn, and 21 



the horn magnetic field, are included in the beam mul- 
tisims. Uncertainties in the magnetic field horn current, 
skin effect of the horn, and secondary nucleon and pion 
interactions in the Be target and Al horn are obtained 
from previous MiniBooNE analyses [ijj ■ 

The normalization of the neutrino component in the 
antineutrino beam was adjusted based on direct mea- 
surements in MiniBooNE [3, [27j ■ The beam fraction of 
Vfi in the antineutrino beam was determined using three 
methods: a pure data sample of events from CCItt 
interactions, differences in Michel electron rates between 
final state [i~ and /i + from and interactions, respec- 



7 



i.600 

j 

;500 

400 
300 
200 
100 







— i — i — i — 


— i — i — i — 


— i — i — i — 


— i — i — i— 


— i — i — i — 


— i — i — i— 


— i — i — _— 


-| IT 


: 
















: 


















: 




























































. . . 






.'~."~T H 









0.4 



0.6 



0.8 



1.2 



1.4 



1.6 



1.8 

= (GeV) 




(a) SciBooNE RS Events 



(b) SciBooNE WS Events 



FIG. 7. Reconstructed antineutrino and neutrino energy (E® E ) distributions for selected RS and WS SciBooNE events on 
different target types (hydrogen or carbon) from SciBooNE MC. Total events are represented by the solid line, events with 
interaction on carbon are represented by the dashed line, and events with interaction on hydrogen are represented by the 
dot-dashed line. 
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FIG. 8. Reconstructed antineutrino and neutrino energy (E® E ) distributions for selected RS and WS SciBooNE events for 
different interaction types (CCQE, CC17T, other) from SciBooNE MC. Total events are represented by the solid line, CCQE 
interaction events are represented by the dashed line, CC17T interaction events are represented by the dot-dashed line, and all 
other interaction (CC multi-7r or NC) events are represented by the short-dashed line. 



tively, due to fi~ capture on carbon, and angular distri- : 
bution differences between final state fi~ and [i + from ■. 
and interactions, respectively. Averaging these three : 
methods, the ir + production in the beam MC was scaled : 
by a factor of 0.78 and given a 12.8% normalization un- : 
certainty. Uncertainties on the production of 7r~ from the 
initial p-Be interaction are calculated using spline fits to : 
data from the HARP experiment [H_J. An updated K+ 
production simulation with reduced uncertainties for the 
initial p-Be interaction is used. This update is based on 
a new Feynman scaling fit [28| to recent SciBooNE mea- 
surements [_-9| . The K° production uncertainties for the 



initial p-Be interaction are from the Sanford-Wang pa- 
rameterization covariance matrix [To| . K~ production 
is estimated using the MARS hadronic interaction pack- 
age [1^ due to the scarcity of production measurements 
in the relevant kinematic regions. K~ production cross 
section uncertainties from the initial p-Be interaction are 
given a conservative 100% normalization uncertainty. 
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and uncertainties of Ma for CC coherent 7r interactions, 
Ma for multi it interactions, Fermi surface momentum 
(pf), and NC axial vector isoscalar contribution (As) 
are identical toprevious MiniBooNE and SciBooNE mea- 
surements 0, [l(| ■ The uncertainties for pion absorption, 
pion inelastic scattering, and pionlcss A decay in the tar- 
get nucleus (±25%, ±30%, and ±100%, respectively) are 
treated in the same way as in a previous measurement 0] , 
however they are treated as uncorrelated between Mini- 
BooNE and SciBooNE (unlike all other cross section un- 
certainties). Both the v and v Ma values and their un- 
certainties for quasielastic interactions on hydrogen are 
based on the latest deuterium measurements [3~i| . 

Additional systematic uncertainties are added to ac- 
count for limitations of the RFG model. Such limita- 
tions include the absence of processes such as meson ex- 
change currents and multi-nucleon knockout events (32r - 
HH. A 10% normalization uncertainty is assigned to both 
v and v CCQE interactions on carbon to cover the differ- 
ence between data and prediction in the MiniBooNE 
CCQE measurement. An additional 40% normalization 
uncertainty is placed on v CCQE interactions on carbon 
to cover the discrepancy between the RFG model pre- 
diction for v and recent nuclear models [3214351 ] . An ad- 
ditional 10% normalization uncertainty is added to non- 
CCQE T> interactions on carbon to account for the limi- 
tations of the RFG model for those type of events. 

The full list of beam and cross section parameters for 
MC simulation and its associated systematic uncertain- 
ties are shown in Table |TTJ 



FIG. 9. Antineutrino and neutrino path lengths for Mini- 
BooNE and SciBooNE events from point of production to 51 
interaction in detector, as predicted by the MC simulation. 

52 
53 
54 

1 B. Cross Section Uncertainties 

55 

56 

2 CCQE cross sections on carbon are calculated assum- 57 

3 ing an RFG model with parameters Ma (axial mass) 5 s 

4 = 1.35 and K (Pauli blocking factor) = 1.007. An ad- 59 

5 ditional correction, as a function of Q 2 , is applied to m 
e background CCItt interaction events in MC [l4j|. The un- 6 i 
7 certainties in Ma and n for CCQE events on carbon arc 62 
a based on the statistical uncertainties of the MiniBooNE 6 3 

9 neutrino mode measurement fl4j , to avoid double count- 64 

10 ing systematic uncertainties accounted for in this analysis 6 5 

11 as detailed in this section. 6 6 

12 Since the purpose of the Q 2 correction in the Mini- 67 

13 B00NE neutrino mode measurement [l4[ is to match 68 

14 the background CC17T interaction events in MC to a se- 69 

15 lected data sample comprising mainly of CC17T interac- 70 

16 tion events, there is no uncertainty placed on Ma for 71 

17 CC17T interaction events. However, for CC17T inter- 72 
is action events, the M^-resonant and coherent l7r uncer- 73 

19 tainties are not constrained by the MiniBooNE neutrino 74 

20 mode measurement and are not reduced. The values 75 



C. Detector Uncertainties 

Uncertainties associated with the MiniBooNE detector 
include light propagation, attenuation, and scattering in 
the detector as well as PMT response. The optical model 
for light propagation in the detector [22j uses 35 param- 
eters for properties such as refractive index, attenuation 
length, scintillation strength, etc. These parameters are 
tuned to non-MiniBooNE measurements as well as Mini- 
BooNE internal data. Over 100 separate MC data sets 
were created based on variations in these parameters. In 
a manner similar to the multisim method, these results 
were used to compute the optical model error matrix in 
bins of reconstructed antineutrino energy. To estimate 
the impact of uncertainties in PMT response, indepen- 
dent MC data sets based on variations in the discrimina- 
tor threshold, or the PMT charge-time correlations, were 
created and compared to default MC. Based on com- 
parisons with external data [36|-[38[ and the output of 
GCALOR, an uncertainty of 35% is assigned to pion ab- 
sorption and 50% is assigned to charge exchange in the 
detector medium. This is distinct from the uncertainty 
on pion absorption and charge exchange inside the nu- 
cleus. 

Uncertainties associated with the SciBooNE detector 
include uncertainties in the muon energy loss in the scin- 
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tillator and iron, light attenuation in the wavelength i 
shifting fibers, and PMT response; see Ref. [l6|. The < 
crosstalk of the MA-PMT was measured to be 3.15% for < 
adjacent channels with an absolute error of 0.4% [12]. < 
The single photoelectron resolution of the MA-PMT is < 
set to 50% in the simulation, and the absolute error is s 
estimated to be ±20%. Birk's constant for quenching in , 
the SciBar scintillator was measured to be 0.0208±0.0023 , 
cm/McV [l^ . The conversion factors for ADC counts to , 
photoelectrons were measured for all 14,336 MA-PMT 
channels in SciBar. The measurement uncertainty was at 
the 20% level. The threshold for hits to be used in SciBar 
track reconstruction is 2.5 photoelectrons; this threshold 
is varied by ±20% to evaluate the systematic error for 
SciBar track reconstruction. The TDC dead time is set 
to 55 ns in the MC simulation, with the error estimated , 
to be ±20 ns [11]. ] 
The reconstruction uncertainties consist of antineu- 
trino energy reconstruction uncertainties and muon track 
misidentification uncertainties. For antineutrino energy 

f 

reconstruction uncertainties, the densities of SciBar, EC, 
and MRD arc varied independently within their mea- 
sured uncertainties of ±3%, ±10%, and ±3%, respec- 
tively. Misidcntified muons stem mainly from proton 
tracks created through NC interactions, which are given 
a conservative ±20% normalization uncertainty. A con- 
servative ±20% normalization uncertainty is applied for 
the MC simulated background of neutrino and antineu- 
trino events initially interacting outside the SciBooNE 
detector that pass the selection criteria. A conserva- 
tive ±20% normalization uncertainty is applied for the 
MC simulated background of neutrino and antineutrino 1 
events initially interacting in the EC /MRD detector that i 
pass the selection criteria. -, 



D. Error Matrix 

All of the MiniBooNE uncertainties, the SciBooNE un- ■ 
certainties, and the correlations between them are ex- ; 
pressed in the total error matrix, M, a 42 x 42 covari- 
ance matrix in MiniBooNE and SciBooNE reconstructed 
antineutrino energy bins defined as: 1 

^MB-SB j^SB-SB 
jyMB-MB j\^SB-MB 



M 



(2) 



4i where 



M 



jur-X. N Y RS N Z RS 

M i,j;(RS,RS) Iy i 



M, 
M 



x 



i.j;{WS,WS) ly i 



jyY WSjqZ WS 



X 



i,j;(RS.WS) ly i 



n7 rs n. 



3 

Z WS 



jT f X N Y WS N Z RS 

M i,j;(WS,RS)^i 

71 ^"X stat 



(3) 



are the bin to bin covariance elements of the full error ma- 
trix. X denotes the type of correlation with Y and Z de- 
noting the type of bins (either MiniBooNE or SciBooNE) 



associated with X. For MiniBooNE to MiniBooNE cor- 
relations, X=MB-MB, Y=MB, Z=MB. For SciBooNE 
to SciBooNE correlations, X=SB-SB, Y=SB, Z=SB. 
For MiniBooNE to SciBooNE correlations, X=MB-SB, 
Y=MB, Z=SB. For SciBooNE to MiniBooNE correla- 
tions, X=SB-MB, Y=SB, Z=MB. N? RS (Nf RS ) and 
n y ws ( N z ws^ arc thc num ber of RS and WS events 
for bin type Y (bin type Z) in reconstructed antineutrino 
energy bin i (bin j), respectively. M^.^ RS RS ^ are the el- 
ements of the RS to RS correlated fractional error matrix 
for correlation type X defined as: 



M i^j;(RS,RS) 



Ivl i,r,(RS,RS) 

jyY RS ]yZ RS 



(4) 



where M^.^ RS is the full RS to RS reconstructed an- 
tineutrino energy bin covariance for correlation type X. 



M x 

1V1 i,j;(WS,WS)> 



M X 

1V1 i,j;{RS,WS)i 



and M lrXWS..RS) 8 



larly defined fractional error matrices for correlation type 
X with different RS and WS correlations. M x stat is thc 
statistical covariance matrix in reconstructed antineu- 
trino energy bins for correlation type X (only SB-SB and 
MB-MB have nonzero elements). 

The decomposition and reconstruction of the full error 
matrix M to and from the fractional error matrices allows 
the error matrix to be updated based on different MC 
predictions, as a function of the oscillation parameters in 
the physics parameter space. 

For illustrative purposes, Figure [TU] shows the square 
roots of the elements of the total fractional error matrix, 

are the elements of 



'M i:j = y/Wj/y/W^i where M t , 
the total error matrix and Ni (Nj) is the MC prediction 
for reconstructed antineutrino energy bin i (j). Figure [TT1 
shows the correlation coefficients of the total error matrix 
in reconstructed antineutrino energy bins. 

Figure[12]shows the MiniBooNE and SciBooNE default 
MC E® E predictions for RS and WS events with error 
bars corresponding to the ^/Mu values of the error matrix 
diagonal elements. 



VI. ANALYSIS METHODOLOGY 

Oscillation predictions are based on a two-antincutrino 
oscillation model, where the oscillation survival probabil- 
ity for a in the beam is given by 



sin 2 26>sin 2 . (5) 

E 

L is the path length in kilometers, E is the antineutrino 
energy in GeV, 9 is the mixing angle, and Am 2 is the 
difference in the squares of the masses of two different 
mass cigenstates. 

The x 2 statistic is formed 



£ (A — Ni) (M _1 ) . . (Dj — Nj), 



(6) 
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TABLE II. Summary of beam and cross section parameters for MC simulation with its associated systematic uncertainties. 



Beam 



Uncertainty 



7r + production in antineutrino beam (from WS neutrino background) 

7r~ production from p-Be interaction 

K + production from p-Be interaction 

K° production from p-Be interaction 

K ~ production from p-Be interaction 

Nucleon and pion interaction in Bc/Al 

Horn current 

Horn skin effect 



12.8% normalization uncertainty [15] 
Spline fit to HARP data 
Table IX in Ref. [28] 
Table IX in Ref. [10] 
100% normalization uncertainty 
Table XIII in Ref. [10] 
±1 kA 

Horn skin depth, ±1.4 mm 



Cross Sections 



Uncertainty 



CCQE Ma on carbon target 

K 

CCQE Ma on hydrogen target 

CC resonant n Ma 

CC coherent n Ma 

CC multi 7r M A 

Eb 

PF 

As 

CCQE on carbon 

CCQE on carbon {v^) only 

non-CCQE on carbon only 

7r absorption in nucleus 

7r inelastic scattering 

7r-less A decay 



1.35 ±0.07 GeV 
1.007 ±0.005 
1.014 ±0.014 GeV 
1.1 ±0.275 GeV^, 
1.03 ±0.275 GeV"3. 
1.3 ±0.52 GeV ~ 

±9 MeV 
220 ±30 MeV/c 

0.0 ±0.1 
±10% norm error 
±40% norm error 
±10% norm error 
25% 
30% 
100% 



This uncertainty is not applied to CCl7r events that are Q 2 corrected. 




5 10 15 20 25 30 35 40 



MiniBooNE SciBooNE 




MiniBooNE SciBooNE 



FIG. 10. Bin-wise square root of the total (statistical and 



systematic errors combined) fractional error matrix y Mij = 



yj Mij I ^NiNj, where Mij is the total error matrix and Ni 
(Nj) is the MC prediction for reconstructed antineutrino en- 
ergy bin i (j). Bins 1 through 21 are MiniBooNE, bins 22 
through 42 are SciBooNE. 



FIG. 11. Correlation coefficients of the total (sta- 
tistical and systematic errors combined) error matrix 
(pij = Mij / (au<Jjj)). Bins 1 through 21 are MiniBooNE, bins 
22 through 42 are SciBooNE. No bins are anti-correlated. 



where (M is the ij'-th element of the inverse of the 
error matrix M, the covariance matrix in MiniBooNE 
and SciBooNE E$ E bins described in Eq.H D { (Dj) is 
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FIG. 12. RS (solid line) and WS (dashed line) contributions to the default MC prediction for MiniBooNE and SciBooNE 
reconstructed antineutrino and neutrino energy (E® E ) distributions. Error bars are the systematic uncertainties from the 
diagonals of the error matrix (an = y/Mii) and do not account for correlations. 



the data count in bin i (j) and Ni (Nj) is the MC pre- : 
diction for bin i (j), in MiniBooNE and SciBooNE E® E 
bins. Ni is the sum of neutrino (WS) and antineutrino 
(RS) events in the i-th bin: 



N t = Np s (Am 2 , sin 2 2 



e)+N, 



(7) 



As shown in Eq. [JJ only the predicted RS event rate : 
depends on the oscillation parameters, Am 2 and sin 2 29, ■ 
for this two-antineutrino oscillation model. The WS flux 
is assumed to not oscillate. The index runs from 1 to 42 
(21 MiniBooNE E$ E bins and 21 SciBooNE E$ E bins). . 
For the physics analysis fitting, a A\ 2 test statistic is 
used as defined by 



Ax 2 = X 2 (N(6 phys ),M(9 phys )) 
-X 2 (N(6 BF ),M(6 BF )) 



(8) 



where 6*bf refers to the oscillation parameters at the best 
fit point and # p hy S refers to the oscillation parameters at « 
a given test point. « 

The method of Fcldman and Cousins is used to 43 
determine the A% 2 value at each point that corresponds 44 
to a certain confidence level of acceptance or rejection. 45 
To obtain the 90% confidence level exclusion region for 46 
disappearance, a A\ 2 distribution is formed for each 47 
point 6>ph ys in parameter space using many iterations of 4s 
generated fake data at that ^phys- The A^ 2 value from 49 
actual data at each # p hy S is then compared to the fake 50 
data A\ 2 distribution at each 6* p hy S - If the A^ 2 value 51 
from actual data is larger than 90% of the all the fake 52 
data Ax 2 values at a # p hy S point, then the 6> p h ys point in 53 
parameter space is excluded at 90% confidence level. The 54 
aggregation of all the excluded 90% confidence level # p h ys 55 
points forms the 90% confidence level exclusion region. 56 

The full error matrix is used to create the fake data 57 
for the Fcldman and Cousins tests. First, a Cholcsky ss 



decomposition is performed on the error matrix M: 
M = LL\ 



(9) 



where L is a lower triangular matrix and L* is the conju- 
gate transpose of L. Then, a vector u is created, where 
each of the n elements, 42 in total, of u are drawn from 
a Gaussian distribution with mean and variance 1. A 
fluctuated fake data histogram is given by 



A fakc = N d 

e fault 



(10) 



where N^efauit is the default monte carlo prediction as- 
suming an oscillation signal with oscillation parameters 
at point # p hy S - 



VII. RESULTS 

Figure [JjJ] shows the observed event distributions, in 
reconstructed antineutrino energy, for MiniBooNE and 
SciBooNE. The systematic uncertainty shown for the MC 
predictions was computed as just the square roots of the 
diagonals of the total error matrix without correlations. 

Tabic irni lists the event counts in each bin, for data and 
for MC predictions. The listed uncertainties are based on 
the square roots of the diagonals of the total error matrix 
without correlations. (The reported SciBooNE data has 
fractional counts due to the manner in which the cosmic 
ray background is subtracted.) 

A MiniBooNE-only disappearance analysis is included 
to give a sense of what the sensitivity would be with- 
out the inclusion of SciBooNE data. Figure HH shows 
the 90% CL exclusion region and best fit point for the 
MiniBooNE-only disappearance analysis, completed 
using the same methodology as the joint disappearance 
analysis except with the exclusion of SciBooNE data, 
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(a) MiniBooNE (b) SciBooNE 

FIG. 13. Reconstructed antineutrino energy (E® E ) distribution for data events, compared to Monte Carlo predictions, for 
MiniBooNE and SciBooNE. Vertical error bars on data are statistical uncertainty. Shaded error band around simulation is the 
systematic uncertainty computed as the square roots of the diagonals of the total error matrix. 



TABLE III. Observed event counts for each MiniBooNE and SciBooNE data bin, MC predictions, and uncertainty. 



Bin Range 


MB Data 


MB MC 


MB 


SB Data^ 


SB Cosmic 


SB MC 


SB 


(MeV) 






Error ± 




Data 




Error ± 


300-400 


565 


413.5 


111.0 


1077.0 


21.0 


997.3 


136.8 


400-467 


2577 


2139.2 


464.8 


966.8 


89.2 


915.6 


141.9 


467-533 


4433 


4039.9 


802.2 


872.8 


85.2 


834.4 


132.1 


533-600 


5849 


5211.0 


1005.7 


854.4 


72.6 


809.4 


132.2 


600-667 


6411 


5725.6 


1108.7 


856.8 


59.2 


790.6 


137.3 


667-733 


6445 


5778.3 


1130.3 


915.0 


51.0 


781.9 


144.3 


733-800 


6090 


5586.8 


1096.9 


849.8 


52.2 


757.3 


139.5 


800-867 


5678 


5268.3 


1044.8 


876.6 


43.4 


717.1 


138.8 


867-933 


5314 


4826.2 


951.8 


787.0 


39.0 


655.8 


138.0 


933-1000 


4624 


4319.6 


865.1 


688.0 


35.0 


639.7 


129.6 


1000-1067 


4015 


3720.3 


747.2 


628.0 


29.0 


580.2 


125.4 


1067-1133 


3349 


3163.6 


642.1 


569.6 


28.4 


488.7 


105.8 


1133-1200 


2965 


2655.9 


554.3 


496.6 


21.4 


403.9 


92.2 


1200-1267 


2464 


2147.2 


453.0 


377.0 


23.0 


308.4 


74.6 


1267-1333 


1937 


1726.4 


367.8 


273.6 


22.4 


228.4 


53.6 


1333-1400 


1534 


1372.0 


297.9 


178.6 


18.4 


150.0 


37.2 


1400-1467 


1227 


1073.3 


238.1 


111.2 


18.8 


89.4 


23.9 


1467-1533 


859 


820.5 


187.7 


65.4 


17.6 


57.1 


15.0 


1533-1600 


679 


607.2 


145.8 


39.0 


17.0 


33.1 


10.4 


1600-1700 


684 


607.2 


149.1 


40.8 


28.2 


27.6 


9.9 


1700-1900 


610 


560.1 


144.5 


37.6 


39.4 


24.8 


7.8 



b The SB data has its SB cosmic data background removed. 



SciBooNE MC prediction, and SciBooNE error matrix u 
uncertainties in the x 2 statistic. The best fit point is u 
Am 2 = 5.9 eV 2 , sin 2 29 = 0.076. At the best fit point, « 
X 2 = 25.7 (probability 12.4%). For the null oscilla- u 
tion hypothesis, X 2 = 28.3 (probability 13.7%). With is 
A% 2 = 2.6, the null oscillation hypothesis is excluded at ie 
52.4% CL. 

Figure[T5]shows the disappearance limit for the joint 18 
disappearance analysis. For Am 2 = 1 eV 2 and Am 2 = 10 19 
eV 2 , the 90% CL limit for sin 2 29 are at 0.12 and 0.03, 20 



respectively. At sin 2 26 = 1, the 90% CL limit for Am 2 
is 0.17 eV 2 . The best fit point from the joint analy- 
sis is Ato 2 = 5.9 eV 2 , sin 2 29 = 0.086. At the best fit 
point, x 2 = 40.0 (probability 47.1%). For the null oscil- 
lation hypothesis, x 2 = 43.5 (probability 41.2%). With 
Ax 2 = 3.5, the null oscillation hypothesis is excluded at 
81.9% CL. All probabilities in both the MiniBooNE-only 
and joint disappearance analyses arc based on fake data 
studies. 

Figure [TBI shows the data to MC ratios for MiniBooNE 
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FIG. 14. 90% CL exclusion region (solid line) and best 
fit point for the MiniBooNE-only disappearance analysis. 
Also shown is the 90% CL result from the 2009 MiniBooNE 
disappearance analysis [H (dashed line) and CCFR [4l[ (dot- 
dashed line). The expected 90% CL sensitivity band from 
fake data studies is also shown (shaded region); lcr (68%) of 
fake data tests, where the fake data had statistical and sys- 
tematic fluctuations but no oscillation signal, had 90% CL 
limit curves in this shaded region. 



and SciBooNE, as well as how the best fit signal modifies 
the MC predictions. From these ratio plots, it can be 
seen how the best fit signal improves the shape agreement 
between data and MC. Figure [TT1 shows the double ratio 



MiniBooNE data 
MiniBooNE default MC 



SciBooNE data ' 
SciBooNE default MC ' 



fill 



In Figure [T7l any common normalization difference is re- 
moved and the expected result is a value of one. The 
double ratio result agrees well with the expectation ex- 
cept where statistics are small. 



VIII. CONCLUSIONS 

An improved search for disappearance using a two- 
detector combined MiniBooNE/SciBooNE analysis has 
been performed. Previous flux and cross section mea- 
surements, as well as an increased data set, have enabled 
a substantial improvement in the sensitivity to dis- 
appearance. The results are consistent with no short 
baseline disappearance of and we have dramatically 
improved on the excluded regions of the oscillation pa- 
rameter space. MiniBooNE and SciBooNE have pushed 20 
the limit on short baseline disappearance of down to 21 



(a) Linear scale on x-axis. 



1 2 P~r-| 1 1 ^j.-.u.j -1 n — "3 — r 1 r-rrrrr. 




sirf(20) 



(b) Log scale on x-axis. 

FIG. 15. 90% CL exclusion region (solid line) and best 
fit point for the joint MiniBooNE, SciBooNE disappear- 
ance analysis. Also shown is the 90% CL result from the 
2009 MiniBooNE disappearance analysis @| (dashed line) and 
CCFR [HI (dot-dashed line). The expected 90% CL sensitiv- 
ity band from fake data studies is also shown (shaded region); 
la (68%) of fake data tests, where the fake data had statisti- 
cal and systematic fluctuations but no oscillation signal, had 
90% CL limit curves in this shaded region. 



roughly 10%, the region of interest for sterile neutrino 
models. 
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FIG. 16. The ratio, with error bars, between data and default MC as a function of reconstructed antineutrino energy (E® E ). 
The ratio of best fit signal MC to default MC is also shown (dashed line). The best fit results from the joint analysis were used 
to generate the signal MC. The shaded regions are the la band from fake data with statistical and systematic fluctuations but 
no oscillation signal. 
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FIG. 17. The double ratio (Eq. lll[) . with error bars, as a func- 15 
tion of reconstructed antineutrino energy (E® E ). Some of the « 
MiniBooNE and SciBooNE uncertainties cancel in this double 17 
ratio. The double ratio where the MiniBooNE and SciBooNE 18 
signal MC based on the best fit results from the joint analysis 19 
are used in placed of data is also shown (dashed line). The 20 
shaded region is the la band from fake data with statistical 
and systematic fluctuations but no oscillation signal. 
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